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ABSTRACT 



The enrichment of Li in the Universe is still unexplained, presenting various 
puzzles to astrophysics. One open issue is that of obtaining reliable estimates 
for the rate of e~-captures on 7 Be, for T and p conditions different from the 
solar ones. This is of crucial importance to model the Galactic nucleosynthesis 
of Li. In this framework, we present here a new theoretical method for calcu- 
lating the e~ -capture rate in conditions typical of evolved stars. Furthermore, 
we show how our approach compares with state-of-the-art techniques for solar 
conditions, where various estimates are available. Our computations include: i) 
"traditional" calculations of the electronic density at the nucleus, to which the 
e~-capture rate for 7 Be is proportional, for different theoretical approaches in- 
cluding the Thomas-Fermi, Poisson-Boltzmann and Debye-Hiickel (DH) models 
of screening, ii) a new computation, based on a formalism that goes beyond the 
previous ones, adopting a mean-field "adiabatic" approximation to the scattering 
process. The results obtained with the new approach as well as with the tradi- 
tional ones and their differences are discussed in some detail, starting from solar 
conditions, where our approach and the DH model essentially converge to the 
same solution. We then analyze the applicability of both our method and the 
DH model to a rather broad range of T and p values, embracing those typical 
of red giant stars, where both bound and continuum states contribute to the 
capture. We find that, over a wide region of the parameter space explored, the 
DH approximation does not really stand, so that the more general method we 
suggest should be preferred. As a first application, we briefly reanalyze the 7 Li 
abundances in RGB and AGB stars of the Galactic Disk in the light of a revision 
in the Be-decay only; we however underline that the changes we find in the elec- 
tron density at the nucleus would induce effects also on the electron screening 
(for p-captures on Li itself, as well as for other nuclei) so that our new approach 
might have rather wide astrophysical consequences. 

Subject headings: Poisson-Boltzmann, Debye-Hiickel, Thomas-Fermi models - 
first-principles methods - Stars: RGB, AGB - abundances - Nuclear reactions - 
Nucleosynthesis - Electron screening 
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Introduction 



Serious problems affect our understanding of the 



Nucleosynthesis (BBN: see e.g. 



Coc & Vangion 



2010 



i evolution in the Galaxy. Big Bang 
and reference s therein) predicts a 



Li abundance higher than observed in extremely m etal-poor object s ( 



and in low-metallicity Main Sequence (MS) stars (ISpite and Spite 



Bonifacio et al. 



20121 ) 



19821 ). On the other 



hand, in the prese nt interstellar medium 



expected by BBN (ICasuso and Beckmann 



much 7 Li ([Fields et al. 



1994 



Alibes et al 



SM) the 7 Li abundance is higher even than that 



2003). As Galactic Cosmic Rays do not produce 



20021 ) we should rely on stellar nucleosynthesis for 



explaining this increase. 

7 Li is however a very fragile nucleus, so that its photospheric concentration in stars is 
easily destroyed when convective processes can carry it to moderately high temper atures (a 



1986al ). Special 



few millions K), where it undergoes proton captures (IBoesgaard fc Tripiccol 
conditions are therefore required for Li production, which are limited to a few astrophysical 
scenarios. These mainly include novae and intermediate mass stars (IMS) undergoing H 
burning at the base of their envelop e through the so-called Hot Bottom Burning process 
(see e.g. 



"ir 



T 



D'Antona &: Ventura! 120101 ). In current stellar models for low mass stars (LMS: 



1997 



; destroyed a 


ready 


Sestito et al. 


2006) 



those with masses below about 2 — 3 M ) Li is ins tead predictec 
in the early phases of evolution, preceding the MS ( iPinsonneaultl 
In stars of mass lower than solar, whose convective envelopes remain large even during the 
MS, Li destruction continues through this stage; for higher masses, instead, it is predicted 
that external convection shrinks and does not include any more zones hot enough to affect 
Li. Contrary to these model expectations, observations of the Sun and of solar-like stars 
reveal that they undergo extensive Li-depleting processes during central hydrogen burning. 
One of the conse quences is that the solar photosphere is about 100 times less Li-rich 



than meteorites ( lAsplund et al. 



20091 ). Li-destruction processes are known to occur also 
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at intermediate effective temperatures, in MS stars o 



so-ca lled Li-dip (see e.g. 



1998|). 



Boesgaard fe Tripicco 



1986b 



the Galactic disc, generating the 



Balachandran 



1995 



Boesgaard et al 



The above phenomena, unexpected from canonical stellar models, have been interpreted 



(sometimes only qua 



convection (IMichaud 



itatiyely) in terms of mixing episode s of a nature different from pure 



1986 



Michaud fe Charbonnau 



induced ef fects (|Charbonnel fe Lagarde 



proce sses (lEggenberger et al. 



199ll ). often attributed to rotationally- 



20101 ) . atomic diffusion, or mag netic dynamo 



2010), like those observed in active stars ( 1 Andrews et al. 



19881 ). However, the very large Li destruction in the Sun and in other main sequence stars 
is still not accounted for in detail, being too small in the models. In more advanced stages, 
i.e. along the Red Giant Branch (hereafter RGB) an d along the Asymptoti c Giant Branch 



(here after AGB) similar phenomena must be acti ve ( 



20121). perhaps again re lated to magnetic effects (IBusso et al. 



Denissenkov et al. 



Palmerini et al 



2007 



2011a 



Majorca et al. 



Nordhaus et al 



20081 : 



20091 ). The consequence is that, in most evolved stars, Li is further 



Palmerini et al. 



2011bl . and references quoted 



depleted as compared to MS stages (see e.g. 
therein) . 

At odds with our need for finding sites where Li is produced in st ars, only few 



red giants (~ 2%) show Li enhan c ement at the i r photosphere (see e.g. 



Charbonnel fe Balachandran 



Lebzelter et al. 



2000 



Kumar et al 



2011 



Brown 



Uttenthaler fe Lebzelter 



1987; 



2010 



20121 ). Their abundances might in principle be produced by coupled 



mixing and nucleosynthesis episodes, in which the depletion of Li by downward diffusion 
is over-compensated by an upward transport of 7 Be from burning regions, at a fast 
enough pace that it survives destruction by p and e~ captures, reaches the stellar 



envelo pe and finally decays there, reproducing Li. 



clear (I Charbonnel fe Lagarde 



2010; 



Palmerini et al 



he situ ation is however far from being 



2011bf ). Quantitative modeling is in 
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particular hampered by a poor knowledge of how the rate of Be decay changes in the rapidly 



varying conditions below the enve 



from the works done for the Sun (jBahcall 



opes of red giants. In fact, extrapolations of this rate 



1962 



IbenetaL 



19671 : 



Bahcall & Moeller 



19691 ) 



are extremely insecure, due the ambient conditions of H burning in evolved stars, which are 
very different from solar. Indeed, in the layers above the H-burning shell of a red giant the 
temperature (T) spans a range from 70-80 MK down to a few MK, while the density (p) is 
lower than in the solar center from one to five orders of magnitudes. 

In general, the pioneering works of the sixties for the Sun were performed considering 
the ionization degree of Be through the Saha equation, and including the contribution 
of free electrons on the assumption that inside a Debye radius around a Be nucleus 
they behave as a Maxwellian gas, thus following a t reatment of the Coulo mb screening 



in the plasma originally due to Debye and H iickel (Debve &: Huckel 



more recent and ge neral approaches do exist (IGruzinov &: Bahcall 



1997 



Sawyer 



1997 



1923) . Although 



Brown & Sawyer 



20111 ) s ometimes demonstra ting the limits of the Debye-Htickel (hereafter 



DH) approximation (j Johnson et al.lll992l ). they maintain a treatment similar to the 
Born-Oppenhemier approach, where the electronic response is considered to be much faster 
than the ionic one. These works almost invariably find for the Sun results very similar to 
those of the classical studies cited above, but they do not consider the physical conditions 
prevailing in evolved stars. 

It has to be further noted that, at the high temperatures of shell H-burning in red 
giants, recombination might occur in states that are highly excited and very close to each 
other. In such conditions, the Born-Oppenheimer approach is questionable; moreover, the 
conditions for the classical DH approximation often do not actually hold (see next sections) 
and one does not know whether this introduces small or large deviations in the capture rate. 



We decided therefore to explore the problem of e captures on 7 Be for the typical 
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T and p values of H-burning layers in evolved stars, both following the traditional DH 
approach and introducing a new treatment, in which the assumption that electrons follow 
a Maxwell-Boltzmann energy distribution is relaxed (considering them, more generally, as 
a Fermi gas) and in which a mean-field adiabatic approximation to the scattering process 
is used. The primary scope is to provide the missing weak-interaction input data for Li 
nucleosynthesis calculations, clarifying whether such data can be deduced in the traditional 
way or not. Subsequently, we also plan to apply our results to a re-evaluation of the problem 
of electron screening in stellar plasmas (see some comments on that issue in sections 3.1 
and 5). 

We remind that other cases exist, involving nuclei heavier than Li, in which the 
predictions of nucleosynthesis are still unreliable for the lack of knowledge about the 
dependence of electron captures on the ambient conditions, during the transport of newly 
produced nuclei to the stellar surface. Crucial isotopes subject to these uncertainties are, 
e.g., 41 Ca and 205 Pb, which are important clocks for da ting the latest nucleosynthesis 



processes before the contraction of the Solar Nebula (see iBusso et al. 



2003 



Busso 



2011 



and 



references therein). If one considers also ft decays, then it turns out that several reaction 
branchings along the s-process path are still affected by our poor knowledge of weak 
interactions in stars, in contrast with the high accuracy of the competing neutron-capture 
processes. 

In general, a better knowledge of radioactive decays would be relevant for a large 
number of physical problems, well beyond the borders of stellar astrophysics. Accurate 
decay rates are e.g. of paramount importance in Earth and Planetary sciences, in order to 
date geological and astronomical processes by estimating the amount of a given long-lived 
species remaining in a sample (a rock or a stellar photosphere). Furthermore, nuclear 
decay provides an impressive source of heat in any planetary body, including the Earth. 
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It is believed that as much as half the heat measured at the Earth surface, corresponding 
to approximately 21 TW, be due to radioactive proce sses involving 40 K, 232 Th, 235 U 



and 238 U, occurring both in t he crust and in the core (iPollack et al 



Lee Sz Steinle-Neumann 



1993 



Fowler 



2005 



20081 ). Heating from shorter-lived radioactivities like Al is then 



believed to provi de the 



the Sun (see e.g. 



Busso 



energ y for melting and differentiating the early solid bodies around 



2011 



and references therein). In a forthcoming paper our results 
will therefore be extended to other nuclei and applied to the clarification of a few such 
problems. 

This paper is structured in the following way. In section 2 we provide a general 
introduction to the current problems of the electron captures on 7 Be. In section 3 we 
discuss the theoretical framework used for solving the ( 7 Be + e~) scattering event, referring 
to Appendix A for a more detailed technical discussion. In the same section 3 we present 
calculations of the electron density at the 7 Be nucleus, obtained by applying several models 
available in the literature (as detailed in Appendix B) and we show comparisons among 
such values. Section 4 then illustrates in a preliminary way some consequences for the 
problem of Li production and destruction in evolved stars, while tentative conclusions are 
drawn in section 5. 



2. The electron captures on 7 Be: state-of-the-art. 



The driving force responsible for the electron capture decay is the weak nuclear 
interaction, a process of very short range that was therefore long believed to be insensitive 
to extra-nuclear factors, notably the chemical environment, the ionization degree, the 
pressure and temperature conditions. Contrary to this simple yiew, many authors report ed 



1972 



Hahn et al. 



evide nce of cha nges in nuclear decay rates wit h temperature (lEmery 
1976 ). pressure ( Lee fc Steinle-Neumannl 2008 ). and the chemical environment (Ra.v et al 
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2002 



Qhtsuki et al 



2007 



Morisato et al. 



20081 ). which are believed to be connected to 



the modification of the electron density at the nucleus, p e (0), induced by the change of 
these parameters. The rate indeed depends on the s-type atomic-orbital wavef unctions, 
the only ones with finite density at the nucleus. Therefore, factors affecting this density 
can appreciably modify the decay rate. Besides T and p values, also the pressure, the level 
of ionization, and the presence of other charged particles, screening the interaction, might 
in principle be of relevance. On the other side, ab-initio calculations of the radioactive 
decay rate at room temperatu re for 7 Be (and several other isotopes) performed by 



Lee fc Steinle-Neumannl (120081 ) showed a very small dependence (within ~ 0.1 — 0.2%) on 
the chemical environment and on the pressure, up to 25 GPa. 

Furthermore, very little is still known about decay events that occur at very high T in 
ionized media for changing p, such as those found in stellar interiors, and similar studies 

theory and intensively debated. T he recent recomme n datio ns by 



Adelberger et al. 


Ib en et al. 


( 


1967) 



(120 111 ) are mainly based on the work by 



Bahcall fc Moellerl (119691 ) and 



(119671 ) . In these works a partial ionization of 7 Be in the Sun was assumed, 



thus the rate now currently used includes contributions from both bound and continuum 



states and the tota 
by 



to-con tinuum capture ratio is 1.217. Different results were presented 
Shaviv fc Shavivl (120031 ). who assumed that 7 Be is fully ionized in the solar plasma; this 
implied an increase of the 7 Be lifetime by ~ 20 to 30% as co mpared to previous r e comm en- 
dations. An even more complicated situation is depicted by iQuarati fc Scarfond (120091 ). 
who recently found a 7 Be lifetime shorter by about 10%, using a modified DH screening 
potential. The present situation is therefore quite unsatisfactory and the uncertainties 
affecting Li abundances in stars have to cope also with this poor understanding of the basic 
nuclear input data. 



This is actually the main motivation of our present attempt, as essentially no one of 
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the existing reaction rate compilations for stellar phys ics can be safely extrap olated from 



the Sun to other situations. For example, the rate by 



Adelberger et al 



(120 111 ) is derived 



from a simp 



Fowler et al. 



e fit o ver a very small domain of the parameter space, while the formula by 



( 119751 ) imposes a specific choice for the density. 



In this context, our goal is to lay the foundations of a theoretical and computational 
method for studying the electron capture and decay rates at high temperature in a 
density-varying medium, to go beyond the existing treatments. We shall probe its validity 
through a comparison with more classical approaches, and shall use the results to reconsider 
briefly the problem of Li production and destruction in red giant stars. 



3. Dynamics of electron capture from ab-initio calculations: our model 

In this section we describe our first-principle approach for computing electron captures 
over a much wider range of T and p than so far possible. These two parameters alter 
considerably the balance between bound and free electronic states contributing to the 
capture. 

While we believe that our formalism is totally general and can be applied even to other 
systems, as a test case of our method we will estimate the decay rate of 7 Be by electron 
captures: 

7 Be + e~ -> 7 Li + v e (1) 

7 Be can decay into 7 Li through different decay channels. In particular, the decay from the 
ground state of 7 Be can occur to both the ground and the first excited state of 7 Li. Of 
course the phase space and the kinetic energy of the neutrinos will be different for these 
two decay paths. At ambient conditions, i.e. for negligible kinetic energy of the impinging 
electron, 7 Be decays in 53 days into the ground state of 7 Li (3/2~) for the 89.7% of cases, 
while for the remaining 10.3% it decays into the first excited state (l/2~), as discussed by 
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Mathews et al.l (119831 ) . The energy of this latter is 477.4 KeV higher than for the ground 
state. The weak branching ratio (BR) at room temperature is therefore 8.709. By defining 
Qo and Qi as the kinetic energies of the neutrinos, escaping respectively from the 7 Li 
nucleus in its groun d and in its first exci ted state, we have Qo = 861.6 keV and Qi = Qo - 



477.4 = 384.2 keV flFuller fc Smith 



20101 ). Being the kinetic energy higher in the first case, 
the available phase space will be larger. Roughly we can estimate that, at a temperature 
T = 10 7 K, BR= 89.7/10.3 x (Q + kT) 2 /(Q 1 + kT) 2 /(Q 2 /Qj) = 8.684. The percentage 
variation of the BR due to an increase of the temperature by five orders of magnitude is 
thus only 0.3%. Of course, even the weak matrix elements will be different for these two 
channels and this should be in principle taken into account. However, our main goal here 
is to estimate accurately the total electron capture decay rate, as this is what we need 
for assessing the astrophysical consequences of our new results. Thus, in the following 
discussion we shall assume that the decay occurs only to the ground state of 7 Li. 

The framework within which we shall calculate the decay rate is given by the theory 
of scattering under two potentials: V, representing the screened, short-range Coulomb 
potential and W, which represents the weak interaction coupling the Coulomb distorted 
initial state and the final channels. We define 0« iP as a free plane-wave, 4>^ k the perturbed 
'in-state', described by a Coulomb-distorted plus an outgoing spherical wave, <pj k and if)J k 
the perturbed 'out-states', which describe asymptotically the emission of a neutrino with 
relative momentum k, released into the final channel / of the target, respectively with and 
without the weak coupling. Then we can write the cross section of the electron capture 
process as: 

Ei, Ef represent the internal energies of the target ( 7 Be) and of the final decay product 
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( 7 Li), p = m e v and k are the relative electron and neutrino momenta in the initial and final 
channels, and v is the electron velocity in the initial channel, relative to 7 Be. In Eq. ([2]) the 
matrix element (0/ J^|0i, P ) must vanish, as the Coulomb interaction does not couple the 
initial and final decay channels ((4>J k\V\4>i,p) = 0). Thus we can define the T-matrix of the 
weak interaction as 

telW+ fe ) = felTVfe) ( 3 ) 
From Eq. ([2]) one can obtain the electron capture rate by multiplying the cross section for 
the electron current: 



2n Wf lfel^lC>f* + Ei - E f - ck ) (4) 



The calculation of the nuclear T-matrix elements, particularly with the inclusion of the 
p 2 /2m e dependence, is very difficult from first-principles and is out of the scope of this 
paper. Indeed p 2 /2m e depends on the temperature of the system. However, the full 
calculation is not needed in our case, and two approximations will enable us to simplify 
the equations without sacrificing accuracy. Firstly, we will neglect the dependence on the 
temperature of this matrix element, and consider only one initial state. This approximation 
is motivated by the fact that the fir st nuclear excited state of 7 Be is at 429.2 keV above the 



ground state ([Fuller fc Smith 



20101 ). corresponding to a temperature of 5xl0 9 K, larger by 
factors of 50 to 100 with respect to the maximum values found during H-shell burning in 
the RGB or AGB phases. 

The second approximation is to model the weak interaction, owing to its very 
short-range nature, by a Fermi contact interaction T\y oc 5{r), independent of the neutrino 
momentum k. By integrating out in k using the latter assumption: 

r *-/ = ^ o|^(o)>| a = ^ |t A4 | 2 <i, o|0+,> + ^ - <0£,|<, o> 

= -L |^| 2 (z,O|0+ p > (Ho + V + Et-Ef) 2 (0,0) (5) 
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where k = ~ ^2m~ + Ei — EA and (0^1^,0) is the electron wavefunction representation at 
the Be nucleus. We will assume that the nuclear T- matrix elements tfj are known and 
equal to those measured on the Earth. The assumptions of Eq. (jSJ) imply that the electron 
capture can be modeled as a 7 Be-e~ two-body scattering process at a given relative electron 
momentum p, and that the rate is proportional to the electron density at the nucleus p e (0), 
which is screened and modified by the presence of the surrounding particles. 

Our model system of the hot plasma, for different conditions of T and p, is represented 
by a homogeneous Fermi gas composed by 7 Be atoms surrounded by N p protons (hydrogen 
nuclei) and N e electrons. The presence of other particles, such as helium, will be neglected 
in the following discussion and in the stellar decay rate calculation, unless otherwise stated. 
However, the generalization of our method to include several species is straightforward. 

An exact, ab-initio calculation of the electron capture rate would be extremely 
complex, due to the many-body nature of the scattering, and one needs to introduce 
some further approximations. The first step of our method is thus the formal reduction 
of this complicated many-body problem to a screened two-body scattering problem, by 
using an "adiabatic" factorization of the eigenfunctions, resembling the widely known 
Born-Oppenheimer (BO) approximation adopted in standard electronic structure methods. 
This goal is reached by fixing the reference frame into the 7 Be nucleus and by writing 
the different parts of the Hamiltonian in this non-inertial frame (we remember that the 
7 Be nucleus is in principle free to move around). In every sense this is only a coordinate 
transformation. However, as in classical mechanics, the consequence of using a non-inertial 
frame will bring about some complicacy (apparent forces), e.g. in this case we will have a 
two-body kinetic energy operator. In the appendix we describe the mathematical details of 
this derivation, from which we obtain that the many-body scattering Hamiltonian in the 
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coordinate system relative to the 7 Be nucleus can be written as: 

3=1 v e " e/ J=l v P j=l 1 e J' J=l 1 f 



Be 



N e N p N e N e N p N p 

~ _ /?' I + _ r ' I + 12 Tr^ _ /?' I _ ?m d ^ 7 Be 

j=1J=1 \ r eJ **p,j\ j=lk=j+1 \ r e,j r e,k\ J=1 K=J+1 l«p,J * M Be 

- E ■ v;, - £ (— v;, • v;,, -=^££ ■ r„ 

J,J'=l v ±>e 7 3j'=l \ & e /He j=1 J=l 



(6) 

where, in order of appearance, one has the electron and proton kinetic energies, the 
electron-Be and proton-Be potential energies, the electron-proton, electron-electron, 
proton-proton interaction, the Be kinetic energy and, finally, terms coupling the different 
particle species. In the inter-particle coupling terms, resulting from our coordinate 
transformation, {r' e j, R' p j} identify the coordinates of the j-electron and J-proton relative 
to the 7 Be coordinate, R^ e = R Bc - 

We then look for separable eigensolutions of the form: 

* (R' Be , {r' e } , {R' p }) = ({r' e },{R' p }) (7) 

V'eeX^Be) = k X (R' Be ) (8) 

This wavefunction factorization differs from the usual formulation of the BO approximation 
in two ways. At variance with the BO approximation, the function $ ({r' e } , {-Rp}) in 
Eq. ((7j), written in the 7 Be reference system, does not depend parametrically on the 7 Be 
coordinates and thus needs to be calculated only once. Furthermore, the BO approximation 
is applicable only when the electronic potential energy surfaces are well separated. These 
are not our conditions, as the electrons, for T and p values pertinent to the burning regions 
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of stars, occupy either highly excited states of the 7 Be atom or continuum orbits. Thus, we 
are out of reach of the BO scheme, which cannot be rigorously applied. 

However, it can be shown that the last two coupling terms in Eq. (J6]), which are 
crucial for the application of our mean-field treatment of the many-body interaction, can be 
neglected. In Appendix A we provide a detailed explanation of the conditions where this 
approximation can be used. 

We notice that, by introducing the two-body framework and the relative coordinate 
system, we can achieve two important results. The first one is that in our approach 
the 7 Be nucleus is in principle free to move around (thus overcoming the limitations of 
the Born-Oppenheimer scheme), even though in general its motion is strongly limited 
by the presence of other particles; the second one is that, by neglecting the two last 
terms of Eq. ([6]), the two-body electron- 7 Be density-matrix can be factorized as the 
product of two one-body density-matrices and thus it is possible to introduce different 
schemes of approximations to the many-body interaction, including the Hartree-Fock's 
one (hereafter HF). Therefore, the screening brought about by all the interacting fermions 
of the surrounding environment, which modifies the two-body electron- 7 Be scattering and 
thus the electron-capture rate, can be now taken into account, using standard many-body 
techniques. 



3.1. Screening at different levels of accuracy: Thomas-Fermi, 
Poisson— Boltzmann and Debye— Hiickel 



While the importance of the screening for the assessment of the electron-capture rate 
is well understood, the approaches used so far are all based on the DH approximation; thus, 



their reliabi 



lb en et al. 



ityis not a priori guaranteed in every situation. Within this approximation 



(119671 ). for example, realized that, in solar conditions, the electron density at 
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the nucleus is reduced by t he electronic screening for both bound and continuum electrons. 
Gruzinov fc Bahcalll ( 119971 ) further improved this model by integrating the density-matrix 
equation to treat on similar grounds the bound and continuum electrons, and by including 
via a Monte-Carlo approach non-spherical charge fluctuations induced by the small number 
of ions within the radius Ad of a Debye sphere. 

Despite these improvements, the conclusions drawn by these works rely on the 
assumption that the hot plasma can be modeled as a classical non-interacting electron 
gas, screened by using the DH model, in thermodynamic equilibrium with a heat bath 
at absolute temperature T. The range of applicability of this approximation is based on 
both classical and statistical considerations: for the former ones, we need a large number 
of electrons and a high temperature; for the latter ones, we require a smooth change 
of the potential over a characteristic distance (Ad), which is large as compared to the 
thermal De-Broglie wave-length of the electrons (Add)- However, in the solar case, where 
T ~ 16 x 10 6 K,p = 150g/cm 3 , we have, for the electrons, Ad = 0.407, Xdb = 0.352 a.u.. 
The conditions are therefore, already for our Sun, at the limits of validity of the classical 
(Maxwell-Boltzmann) gas approximation. Over the more extended range of parameters 
characterizing the radiative layers above the H-shell in a red giant (down to T ~ 2.0 x lO 6 -^) 
the application of the Maxwell-Boltzmann distribution might no longer be justified. 

We underline that the simple DH approach, commonly used in the literature so far, 
can actually be derived as a two-step approximation of the more general Thomas-Fermi 
(hereafter TF) model, by using the linearized Maxwell-Boltzmann distribution. (In 
Appendix B we provide a detailed discussion of the DH treatment). 

Hence we can compute the crucial parameter, i.e. the electron density at the Be 
nucleus p e (0), at different levels of approximation, in order to disentangle the differences 
introduced in the results by the various approaches. In Table 1 we report the results of 
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such calculations, adopting as an example the solar conditions. We underline that p e {0) is 
directly related to the electron-capture rate. 

In the above calculations, the plasma was modelled by a 7 Be nucleus surrounded by a 
Fermi gas occupying an infinite volume and helium was explicitly included. One can see 
that both the classical Poisson-Boltzmann approach (PB) and the DH approximation give 
sizeable different results (of the order of 38% and 6% respectively) with respect to the TF 
approximation, already for solar conditions. The change becomes even larger when we use 
our more general method, based on the HF approximation, which we will discuss in the 
next section. (Surprisingly, the DH approximation is closer to the HF and TF results than 
the PB case, from which it is derived by linearization). 



Table 1 



Phf(0) 


Ptf(0) 


Ppb(0) 


Pdh(0) 


16.320 ±0.020 


16.135 ±0.025 


14.870 ± 0.020 


16.085 ±0.025 



Table 1: Values (in atomic units) of the electron density at the Be nucleus p e (0) for differ- 
ent levels of the theory (HF=Hartree-Fock, TF=Thomas-Fermi, PB=Poisson-Boltzmann, 
DH=Debye-Huckel) for solar conditions (T = 16MK). The conditions at which we performed 
the calculations are T = 15.67MK, p=152.9 g/cm 3 , X H = 0.34608, X Hc = 0.63368. 

Even more important, one must underline that the PB and DH approximations do 
not hold outside the conditions of the solar nucleus. This is true, in particular, at lower 
temperatures and densities, where a large part of the Li production occurs (because the 
competing p-captures on 7 Be become ineffective). In order to check this point, we performed 
the calculations of the the electron density at the Be nucleus for all the previously discussed 
approaches over a wide range of T and p conditions. The results are reported in Table 
2, where T is in units of 10 6 K and p in g/cm 3 . Values at the left of the vertical bar (|) 
represent the electron density at the nucleus p e (0), while those at the right represent p e (0) 
multiplied by (1 + p 2 /2m) /(Ei — Ef) (see Eq. [5]), to which the capture rate is proportional. 
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The Debye radius Ad and the De-Broglie wavelength \ DB for electrons (e ) and protons 
(p) are given in atomic units. 

From the analysis of the electron density at the nucleus p e (0) for the different 
approximations, it is clear that, when we are out of reach of the DH approximation 
(\ D < Xdb), the values calculated using this model may differ by more than 30% from 
the HF values. In particular, we underline that, except for very high temperatures, the 
PB approximation is extremely unreliable, while over the whole considered range of T and 
p values, the TF model systematically underestimates the density at the nucleus p e (0). 
At variance with this, the DH model, except for low temperatures (T < 10 6 K), seems 
to partially correct the wrong behavior of the PB approximation, but not in a systematic 
way. This is due to the fact that the electron density at the nucleus p e (0) is finite for the 
DH model, while it is infinite for the Boltzmann approximation (in the HF method the 
density is always finite, due to the fact that in this case the electronic wavefunction is 
finite at the nucleus). It is thus evident that the DH model fails in capturing the electron 
screening properties over a portion of the parameter space that is very important for Li 
nucleosynthesis and, if one is to assess accurately the electron-capture rate of 7 Be in that 
region, it is necessary to go beyond this approach. 

We notice that the issue at stake here is not only a better estimate of weak interactions. 
If the traditional methods underestimate the electron density over a considerable part of 
the conditions typical of H-burning in Red Giants, as it seems to be the case, then this 
implies that they also underestimate the effect of electrons in screening strong interactions. 
In the case we discuss here (Li production and destruction) our next step will therefore be 
a re-evaluation of electron screening for proton captures on 7 Li, which occur down to very 
low temperatures (few 10 6 K). It is well known that the huge destruction of Li in the Sun is 
still unaccounted for quantitatively, being too small in the models, so that a more effective 
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process of proton captures on Li, induced by an increased estimate of the electron density 
and of its screening of the Coulomb barrier, might have important consequences on solar 
physics and on our understanding of nucleosynthesis in general. 

At high temperature and high density (T > 10 7 K, p > 100 g/cm 3 ) the approaches 
discussed so far (with the exception of the purely classic PB treatment) become essentially 
equivalent to one another, as shown in Table 1; in particular, the values of the electron 
density at the Be nucleus produced by our new method and the one from the DH 
approximation differ by less than 1.5% (see Table 1, column 1 and 4). 

3.2. Calculations of an accurate decay rate from first principles 

In order to apply our model Hamiltonian, Eq. fl6]), to the calculation of the capture 
rate over a wider range of T and p values than so far possible, we need to develop new 
tools, suitable to go beyond the DH approximation and capable of treating accurately the 
electronic screening for both bound and continuum electrons. We will then use the new 
estimates for interpreting the observational data. 

In our approach, the Coulomb screening is calculated using a temperature-dependent 
HF method, within the canonical ensemble, by populating both the ground-state and 
excited-state atomic orbitals via a Fermi-Dirac distribution. In our treatment, we will 
neglect the electron-proton pair- formation at low temperature (T < 10 6 K). Our method 
can be outlined as follows, i) We assume that the decay rate is proportional to p e (0) 
calculated by solving self-consistently a system of coupled HF equations for protons and 
electrons in the electrostatic field of a 7 Be nucleus at the origin of the coordinate system, 
ii) The coupling is given by the self-consistent Hartree term, describing the Coulomb 
interaction between the two fermionic species, iii) The chemical potentials of protons and 
electrons are chosen in such a way that the system is neutral far from the 7 Be atom, iv) 
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Within this framework, the scattering of 7 Be with protons and electrons is a two-body 
interaction, with the remaining particles acting as a mean-field. 

The first step in the above list consists in solving self-consistently the equations for the 
electron charge density p and the potential V to calculate the electron-capture rate. The 
HF equations for this problem read: 



Pj 



,aa e T Be a'cr' E T Be \ r : r 



E 



e kT + l 



j,aa e T Be 



n{r)ipj 



Ar') 



(10) 



VI- 



/Z Z i 
ds urr^T^^^ e ^ e Pi',/3 CT »^ e /3<^ e (s, s) (ii) 
1 1 



I ' ' 

In Eqs. P EDim^G) and r ^ (I) 

are the electronic and nuclear spins, while a and /3 
represent all the other quantum numbers. The index j runs over all the fermionic particle 
types. The self-consistent HF potential is used in Eq. (jSJ), along with the external potential, 
V = Vhf + Vext, t° calculate the (static exchange) electron-capture decay rate. 

In order to be complete, the above treatment of the many-body interaction needs 
to include the electron and proton continuum states, as the capture can occur from 
the continuum orbitals and, at high T, the plasma is from partially to totally ionized. 
To include the continuum states in the HF equations, we used the theory of projected 



potentials, develope c 



(jTaioli et al. 



2009a.b 



or t 



re calculation of the electronic emission spectra from solids 



20101 ). Within this theory, while the Coulomb interaction among 



the particles (see Eq. [6]) is projected onto the Hilbert space spanned by a basis set, 



-20 - 




Fig. 1. — The electron-capture half-life in days for 7 Be as a function of p and T. 

giving rise to the discrete part of the spectrum, the kinetic terms are left unprojected. In 
this way, as the eigenfunctions of the kinetic energy operator are plane waves at a given 
energy, one can "recover" the energy continuum. Furthermore, the sums over the discrete 
states n in Eq. ( FlOl) should be thought of as integrals over the continuum. In our case a 
cc — pVDZ Gaussian basis set (GBS), centered on the 7 Be nucleus, has been optimized for 
the calculation of the bound states. We underline that the addition of the continuum states 
to the projected Coulomb interaction can have significant enhancement effects on the decay 
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Table 3 





Energy 


Pet(0) 


Hartree-Fock 


-14.5729914739 


17.68521 


PuU-CI 


-14.6604710493 


17.68060 



Table 2: Energy of the isolated beryllium atom in atomic units and spin-up density at the 
nucleus obtained through the HF and CI calculations. 

rate over a wide portion of the parameter space spanned. This however does not include 



the solar case, w 



Adelberger et al. 



r ere ou r estimates are essentially indistinguishable from those reported by 
(120 111) . We notice in any case that, whe n applicable, our c hange s point to 



Shaviv fc Shavivl (120031 ) for the 



the opposite direction with respect to the suggestions by 
Sun. 

Our two-body scattering framework allows us to go beyond even the HF treatment of 
the screening; therefore, we estimated the importance of correlations by computing p e (0) 
for an isolated 7 Be atom using the Full Configuration Interaction approach (FCI) for the 
previously used cc — pVDZ GBS. We did not find any appreciable difference between the 
p e (Q) estimates calculated by the HF and by the FCI methods, up to the third digit (see 
Table 3), thus justifying our mean-field approach, which neglects dynamical correlations. 

The calculation of the decay rate, performed at different T and p values using our 
approach is reported in Figure [TJ We recall once again that the above treatment of the 
many-body interaction goes beyond the DH approximation previously used to estimate 



the electron-capture decay rate in the Sun (jlben et al. 



1967 



Bahcall 



1962|). The field 



of applicability of the HF screening effects can in fact be extended rigorously, to cover 
the whole range of parameters found in the layers between the H-burning region and the 
envelope base in evolved stars. 
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4. An application to the Li production and destruction in evolved stars 



As an example of application of our new rate estimates to practical nucleosynthesis 
problems in stars, we briefly illustrate here the situation for rather low mass stars 
(M < 2M Q ) of solar metallicity, undergoing the RGB and AGB evolu tionary phases 



i n pres ence of deep-mixing processes. They were recently discussed by 



( ]2011al Jbl) within the framework of reaction rates offered by 



Palmerini et al. 



Adelberger et al. 



(2011). 



We consider here only the impact on the above astrophysical scenario of our "best 
choice" for the reaction rate (our new method), without analyzing the individual behavior 
of all the four estimates presented, as we are preparing a more thorough analysis on that, to 
appear in a forthcoming paper, also dealing with various stellar masses and compositions. 
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Fig. 2. — A comparison between the equilibrium abundances of 7 Be and 7 Li achieved in the 
layers ab ove the H-burning shell, adopting our 7 Be life-time (red line) and the one extrap- 
olated by lAdelberger et al.l ( 1201 ll ) (black line), in a 2 M evolved star of solar metallicity. 
Panel a) refers to typical RGB conditions, panel b) to AGB stages. The matter density is 
also shown (blue line), and is referred to the scale on the right axis. 



Figure [2] (a,b) shows the abundances of 7 Be and of its daughter 7 Li in the layers below 
the convective envelope bottom and above the H-burning shell. They are plotted shortly 
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after the bump of the Luminosity Function on the RGB (panel a) and in between two 
thermal pulses, d uring quiet H-burning, o n the AGB (panel b). The stellar models are the 



same discussed in 



Palmerini et al. 



( 1201 lbl ). The red lines illustrate the situation obtained 



with our new rate (the case previously labell ed HF), while the black line shows the previous 



findings, obtained with extrapolations from 



Adelberger et al. 



(120 111 ), hence based on the 



DH approximation. The right-end limit of the plot, at the position Ar = 0, characterized 
by temperatures below 10 7 K, represents the base of the envelope; the left-end limit is the 
region where the maximum energy is released from H-burning. As is shown, the matter 
density of the layers considered never achieves the high values typical of the solar core: this 
is the most critical reason why the DH approximation is inadequate for evolved stars. In 
our more general approach both higher and lower values of the electron density near the 
7 Be nucleus can be found, depending on the conditions; however, over a major part of the 
region of interest, where the density is not far from one tenth solar and the temperature 
remains high enough (T > 10 7 K) the values of p e (0) we find are higher t han in the DH 



Palmerini et al. 



( 1201 lbl ). and a 



model; the electron captures are therefore faster than in 
lower equilibrium abundance of 7 Be is established. 

Notice that in our present example neither the 3 He+ 4 He rate, nor the 7 Be+p or the 
7 Li+p rate have been modified, so that the variations shown in Figure H] are entirely due 
to the new approach adopted in computing electron captures on 7 Be. As discussed before, 
further effects are actually to be expected on proton captures, as the changes in the electron 
density will affect the screening of the Coulomb barrier for charged-particle interactions. 

As discussed in the Section 1, it is known that interpreting the observed isotopic 
abundances of light and intermediate elements in evolved stars does not require only the use 
of the proper reaction rates, but also needs the assumption that deep mixing phenomena 
occur, both destroying fragile nuclei of the envelope (like Li), by exposing them to high 
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temperatures, and carrying to the surface products of H-burning nucleosynthesis. These 



mixing mechanis ms are generally 



the approach by 



Palmerini et al 



param eterized and depend on two main parameters; in 



( 12011a 



these parameters are the mass circulation rate 



(in units of 10~ 6 M /yr) induced by the transport processes (M 6 ) and the depth in the 
structure they achieve. This last is expressed in terms of the (logarithmic) temperature 
difference between the deepest mixed layers and those where the maximum energy is 
released in the H-shell (A = logT# — log Tp). (Recently, an analysis devised specifically for 
fixing s uch extra-mixing parameters, adopting bright RGB stars as constraints, was carried 



out by 



Abia et al. 



2012 



, confirming the indications by Palmerini et al. on the need of a 



rather shallow and slow mixing on the RGB). 




Fig. 3. — Panel a) The Li destruction induced by extra-mi xing on the RGB. Panel b) The 
same effects for AGB stages (see text for explanations, and iPalmerini et al. fl2011bh for the 
details of the observational points and their references) 



With reference to Panel a) in Figure the points labeled BLF indicate the position in 
temperature of the Bump of the Luminosity Function, for two representative Li abundances, 
inside the wide spread allowed both by observations and models of the previous evolutionary 
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stages (-0.3 < v4(Li)< 1.5); crosses indicate observations. From each one of the two BLF 
conditions chosen, the two black curves show the Li destruction induced by adopting, for 
the extra-mixing parameters, A = 0.22 and M 6 = 0.01 (upper curve) or 0.1 (lower curve). 
The effects of deep mixing are considered up to the moment when the RGB tip is reached. 
Subsequently, the temperature increases at constant values of A(Li) (dashed lines) up to 
the start of core He-burning (in the CLUMP region in the plot), after which Li would 



decrea se again during the early- AGB stages (e-AGB in the plot; see also 



Palmerini et al 



2011bl . for details). The curves showing the evolution of Li are completely insensitive to the 
adoption of the new or of the old rate for e~-captures. Panel b) of Figure [3] shows instead 
that a different situation emerges during the subsequent AGB evolution. Here we started 
the computations from three initial values of the Li abundance, inside the observed spread 
left by RGB stages. Open and solid dots in the Figure represent observations of M-MS-S 
and of C(N) stars, respectively. Asterisks show the high Li abundances of CJ stars. The 
curves refer to stellar models with A = 0.22 and M$ = 0.1, 0.3 1 (cases labelled 1, 2, 



3, respectively). 



Adelberger et al. 



he b lack curves are obtained using the 7 Be + e rate extrapolated by 



20111 ). while the red ones refer to our best choice from the present paper 
(HF). It is clear that large changes emerge (in particular, with the new rate, Li production 
becomes impossible). 

The reasons for the above dichotomy between RGB and AGB stages can be understood 
with reference to the time scales for 7 Be mixing and decay in the two cases. They are shown 
in Figure HI As the figure shows, for the range of mixing rates required to explain the other 
RGB chemical anomalies (panel a) both the new and old choice for the 7 Be decay provide a 
short enough lifetime for 7 Be, to allow its decay to Li before being saved to the envelope. 
In other words, conditions suitable for Be destruction are in any case met, both with the 
older and with the newer choice of the rate. This explains why the new estimate for the 
rate, which has evident effects in the stellar structure shown in panel a) of Fig [2], does not 
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modify the envelope abundances produced by the mixing process. The opposite occurs for 
AGB stages (panel b). Here, the reduction in the lifetime of 7 Be obtained with our new 
rate is strong enough that even with relatively fast mixing processes the destruction of Li 
in the envelope is not compensated, and the highest observed Li abundances (typical of 
CJ stars) cannot be reached. Essentially, with our rate Li behaves, on the AGB, in a way 
similar to what occurs on the RGB: destruction always prevails. 

We notice that, far from creating a new problem, this finding is now in agreement 
with the indications provided by other chemic al anomalies of the very peculiar CJ stars. 



Let's comment on this in some more detail. In 



Palmerini et al. 



( 1201 lbl ) it was surprisingly 



shown that, adopting t he extrapolation to AGB conditions of the 7 Be lifetime from 



Adelberger et al. 



(120111 ). the high Li abundances observed in CJ-type carbon stars could 



be explained in the framework of the evolution of si ngle stellar structures e xperiencing 



Palmerini et al. 



201 lbh . This result 



rather fast extra-mixing (see in particular Fig. 10 in 
was in itself quite strange, as many peculiarities in such objects suggested instead, since a 
long time, that they do not follow the "normal" evolutionary sequence for LMS. This is in 
particular demonstrated by their luminosity, on average lower than for other C stars, by 
their lack of enrichment is s-process elements and by the fact that a remarkable fraction of 
CJ stars have O-rich shells. Among the hypotheses presented in the literature for explaining 
the ano malous CJ evolutionary path, binarity is probably the most commonly-invoked 



scenario ( ILambert et al. 



1990 



a single peculiar object ( iMcClure 



Abia k 



sern 



20001 ) . especially if leading to coalescence into 



19971 ). In this case, obviously, the high Li abundances 
of the peculiar descendants of such a stellar merging are not expected to be explained by 
single-star nucleosynthesis. We actually predict that Li destruction, with our new rate, be 
even more facilitated than what emerges from the first simple example of Figure El In fact 
a further effect is expected by an increase in the electron screening, due to the now higher 
(on average) electron densities near the nuclei, making proton captures on Li more effective. 
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In general, the Li-rich red giants (even outside the peculiar CJ class) need really 
exceptionally fast transport rates to explain the net Li production. Their small number says 
that these phenomena must be very rare, or concentrated in extremely short evolutionary 
stages. 




Fig. 4. — A comparison betwe en our life-time for 7 Be, T7 Be (red line) and the one extrapolated 
from Adelberger et alj ( 201l[ ) (black line) with the time scales of non-convective mixing at 
various rates (labelled by the mixing rate in units of 10~ 6 M /yr) in the radiative layers 
between the H-burning shell and the base of the convective envelope, for a star with 2 M Q 
and solar metallicity. Panel a) presents a typical situation for the RGB phase, while panel 
b) is instead typical of AGB conditions. 



5. Conclusions 

In this note we have presented a revision of the methods normally employed for 
estimating electron captures and, more in general, the electron density near a nucleus in a 
stellar plasma. For this scope we have devised an ab-initio technique, based on a formalism 
that goes beyond the previously-adopted ones. In our method, we firstly reduced the 
complicated many-body problem by a screened two-body scattering model, by using an 
"adiabatic" factorization of the eigenfunctions, resembling the Born-Oppenheimer (BO) 
approximation. We did that by fixing the reference frame into the 7 Be nucleus and then 
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writing the different parts of the Hamiltonian in this frame. The frame itself was considered 
as non-inertial, to account for the possibility that 7 Be be in motion in a rather complex 
way. In this coordinate transformation, the non-inertiality translated into some technical 
complicacy, with the need of considering the ensuing apparent forces. The mathematical 
details of this procedure were presented in an Appendix. In our Hamiltonian we then had 
to include the electron kinetic and the electron-Be potential energies, the proton kinetic and 
proton-Be potential energies, the Be kinetic energy, the electron-proton, electron-electron, 
proton-proton interactions. Furthermore, we also had to consider two-particle kinetic 
terms, which identify the coordinates of the j-electron and J-proton relative to the 7 Be 
coordinate. In particular, the conditions under which the two terms coupling the Be 
and plasma momenta must be considered carefully, as this is crucial for the application 
of our mean-field treatment. We then looked for separable eigensolutions and wrote 
the wavefunction factorization in a way different from the usual formulation of the BO 
approximation. Indeed, at variance with that approach, in our method there is neither 
the need that the electronic potential energy surfaces be widely separated, nor that there 
is a mass scale difference among the particles. By introducing the two-body framework 
and the relative coordinate system, we achieved two main results. The first one is that 
in our approach the 7 Be nucleus is free to move around (thus overcoming the limitations 
of the Born-Oppenheimer scheme); the second one is that the two-body electron- 7 Be 
density-matrix can be factorized as the product of two one-body density-matrices and thus 
it is possible to introduce different schemes of approximations to the many-body interaction, 
including the Hartree-Fock's one. Then, the screening effects due to all the fermions of 
the environment, which modify the electron- 7 Be scattering, were taken into account, using 
standard many-body techniques. 

By comparing our results with those obtained at various levels of simplification of 
the theory we showed that, while the traditional DH approximation yields reliable results 



-29 - 



for the solar core, where both high temperature and high density conditions exist, such 
simplified methods do not always hold in stellar evolution. In particular, in the low density 
and temperature environments characterizing partial H-burning above the H-shell in Red 
Giants, where lithium undergoes important depletion, we showed that our general method 
should be used and that extrapolations from existing e~-capture rates on 7 Be in the Sun 
are not applicable. 

As an example of the astrophysical consequences of our new method we recomputed, 
for a star of M — 2M and solar metallicity, the evolution of the surface Li abundance in 
presence of deep mixing mechanisms. While, for RGB phases, with the new technique we 
re-obtain previously known results, on the AGB (where the density drops more remarkably) 
a stronger Li depletion was obtained and it was shown that high Li abundances (as observed 
e.g. in CJ-type carbon-rich giants) cannot be the result of extra-mixing occurring in single 
stars, but need more intricate evolutionary paths, as already speculated in the past from 
other chemical peculiarities of these objects. 

Apart from the limited example presented in this first contribution, our method yields, 
in the regions characterized by temperatures of one to a few 10 7 K and moderate densities 
(from 1 to 30-40 g/cm 3 ) higher values of the electron density at the nucleus than the DH 
approach. These conditions are not relevant for the solar core, but are important in red 
giant stars and in sub-envelope layers of MS stars. This should have the effect of enhancing, 
in those astrophysical environments, the electron screening, thus favoring thermonuclear 
fusion. Such a possibility may be especially important for light nuclei (like Li, or 3 He), 
burning at relatively low temperature and that are related to crucial cosmological problems, 
like the explanation of the gradual 3 He consumption in the Universe and of the huge Li 
destruction below the surface of main sequence stars, including our Sun. 
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A. Transformation to relative coordinates 

The calculation of the electron-capture decay rate in a many-body system is a difficult 
task, thus we propose to use an "adiabatic" approximation of the Hamiltonian, which 
resembles the widely used Born-Oppenheimer approach. 

In order to derive the equation of motion for the fermions in the coordinate system 
relative to the 7 Be nucleus, we will consider a model system composed by a single 7 Be 
atom surrounded by N p protons and N e electrons. The Hamiltonian of such system in the 
laboratory frame can be written: 

N e / 1 v. N b , . >. N e 



'Be 



2M Be 2m e e 'V 2m p P ' J J \R Be - r eJ \ 

N P N e Np N e N e 



E \R Be Rp,j\ ^J^kej Rp,j\ ^^_ 1 \r e ,j r e ,fc| 

N p N p 

+E E lR R | (ad 

J=1K=J+1 \ K P,J- K P,k\ 

By performing the following transformation to the ( 7 Be-e~ and 7 Be-p + ) relative coordinates 

^Be = R-Be, r ' e ,j = r e,j ~ -Rfie, R' p j = Rp,J ~ Rbc 

Np N e 

v Be = v Be - v;, , - £ v;,, v^. = v e „ v; 5j = v PiJ (A2) 
j=i j=i 
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one obtains: 
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3=1 J=l e J P> J| i=l fc=j+l 1 e J e > fcl J=lA-=J+l' P. J P>^ 

- y (— v r -v «] - y f — v -v ,)-— yyv r -v 



RLr-RLJ 2M Be ' Be 



J,J'=1 ^ ±(e / jj'=l v ae / ae j=l J=1 



(A3) 

Within this model the 7 Be nucleus plays the role of the heavy nucleus in the BO 
approximation, moving more slowly than the light particles (electrons and protons). 

Eq. (1A3j) can be simplified by neglecting all the inter-particle coupling terms, which 
are divided by the Be mass (M Be ). However, the last two terms in Eq. (1A3j) are special 
as they couple the Be nucleus and plasma momenta. Thus, we need to find the conditions 
under which these two coupling terms can be neglected as this is crucial for the application 
of our mean- field treatment. To determine when this is possible, one can notice that the 
two-body density-matrix of the system in relative coordinates at temperature T is given by: 

1 f elk ^ r Q k i }P~ 
Z=zl2j (^e~^e~^5^|$ a , fc >< $ Q)fe | <g> \ Xk >< Xk\ (A4) 

where a identifies the electronic quantum numbers, E r a k are the eigenvalues of the secular 
problem for $, and Z is the canonical partition function. The exponential term in Eq. ( 1A4I) 
kills the integral, unless: 

k ~ ^/2M Be k B T (A5) 

Using Eq. (IA5j) . the two coupling terms in Eq. fl6]) are of the order of ■ Vg^| ~ 

2^k B TK e m e /M Be and • V p ,j\ — 2\/k B TK p m p /M Be , where K e and K p are the 
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electron and proton kinetic energies, respectively. Therefore, if the two conditions 

Ak B Tm e /M Be < K e , Ak B Tm p /M Be « K p (A6) 
hold, then ( — —k ■ VI ) and ( — — k ■ V' .■ ) are negligible. These conditions are generally 

\ Trie. 'J J \ Tflp PiJ J 

satisfied, due the presence of the Be mass in the denominator and to the fact that the scalar 
products {^—^-k ■ Vg^-J and {^—-^-k ■ V^jj contain the cosine of the reciprocal direction 
of the two multiplying vectors, which, on average, is very small. Finally, the last three 
terms of the Hamiltonian ( 1A3I) can be safely neglected as they contain a multiplying factor 
proportional to the inverse of the Be mass. 



B. The Thomas— Fermi and Debye-Hiickel models 

The widely used DH model can be obtained via a two-step approximation starting 
from the TF theory of the electron gas. The TF model is a simplified HF theory, in which 
the electron gas is treated within the local density approximation (LDA), so that a large 
number of electrons is needed in a region where the potential is nearly constant. In this 
theory the density is thus diagonal in the electron coordinates: 

Pj, aa '(r,r') = Pj, aa '{r)6(r - r') (Bl) 

/ \ f dp tpjap' l Pj,a'p /t-)o\ 

J [iix) 6 e e J>p + 1 

and can be obtained by the self-consistent solution of the TF equation: 

- dr^W + V £p ( r )VW - A*WW + Yl KcpPp^Pp = e J,P^W ( B3 ) 

where V^{r) is the electron-nucleus interaction, \i is the chemical potential, a, a' identify 
the electron quantum numbers and 

V 3 %{r, r') = 6(r-r') ^ / dsg ja ^ fa ^> (r - s) pf,pp(s, s)-J29jap„jp>a>(r-r')p j> pp l (r, r') 

(B4) 
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. Here g is the bare Coulomb potential. 

In order to obtain the DH approximation, in the first step one substitutes the 
quantum-mechanical Fermi-Dirac statistics with the classical Boltzmann distribution, 
while, in a second step, the high-temperature (weak coupling) limit is obtained by Taylor 
expanding at the first order the exponential distribution: 

In this way one obtains the DH equation for a neutral plasma as follows: 

V 2 $ = A D $ (B6) 

where 

*w = £ / 

j' 1 1 

In Eq. (IB6j) . the DH length is defined by Ad = (e r eo&BT/ Ylf=i Pj<?J)^ where k B is the 
Boltzmann constant, p° is the mean charge density of the species j, e r and eo are the relative 
and vacuum dielectric constants. In this framework, A^> sets the characteristic length scale 
for the variation of the potential and of the charge concentration. 
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Values of the electron density p e (0) (in atomic units) at the Be nucleus, for different theoretical approaches (see text for explanations) 



